Spin polarization in the Hubbard model with Rashba spin-orbit coupling on a ladder 



o 

(N 



u 



I 

o 
o 



> 

rn 
cn 

o 



Jose A. Riera 
Instituto de Fisica Rosario y Departamento de Fisica, 
Universidad Nacional de Rosario, Rosano, Argentina 
(Dated: March 18, 2013) 

The competition between on-site Coulomb repulsion and Rashba spin-orbit coupling (RSOC) is 
studied on two-leg ladders by numerical techniques. By studying persistent currents in closed rings 
by exact diagonalization, it is found that the contribution to the current due to the RSOC Vso, for 
a fixed value of the Hubbard repulsion U reaches a maximum at intermediate values of Vso- By 
increasing the repulsive Hubbard coupling U, this spin-flipping current is suppressed and eventually 
it becomes opposite to the spin-conserving current. The main result is that the spin accumulation 
deflned as the relative spin polarization between the two legs of the ladder is enhanced by U. Similar 
results for this Hubbard-Rashba model are observed for a completely different setup in which two 
halves of the ladders are connected to a voltage bias and the ensuing time-dependent regime is 
studied by the density matrix-renormalization group technique. It is also interesting a combined 
effect between Vso and U leading to a strong enhancement of antiferromagnetic order which in turn 
may explain the observed behavior of the spin-fiipping current. The implications of this enhancement 
of the spin-Hall effect with electron correlations for spintronic devices is discussed. 

PACS numbers: 71.27.+a, 71.70.Ej, 73.23.-b 



Introduction. Recent trends in the field of spintronics 
[ll-Q exploit the possibility of controlling electron spins 
by purely electrical means without using magnetic mate- 
rials [3| , which is at the heart of the conceptual proposal 
of the Datta-Das spin field-effect transistor [5] . Key to 
this possibility is the Rashba spin-orbit (SO) interaction 
@) 0] appearing due to structure inversion asymmetry of 
materials and which has been shown to be tuned by gate 
voltages in semiconductor heterostructures Q. 

Although the main interest in the study of RSOC has 
been on its ability of generating and controlling spin po- 
larized currents, another consequence of the interplay be- 
tween charge and spin degrees of freedom induced by the 
RSOC is a flow of spins transversal to the flow of charge 
which is known as the spin- Hah effect If transport 

takes place along a planar strip, the spin-Hall effect leads 
to a spin accumulation with opposite projection on both 

This has been experimentally 



edges of the strip [ 

observed [l^ Hence, the spin-Hall effect provides 

a real-space separation between up and down spin elec- 
trons, like two ferromagnets with opposite polarizations, 
thus enhancing the metallic state and its conductance, 
and these features could also be employed in spintronics 
devices. 

Since most devices and experiments about RSOC 
have involved semiconductors, electron correlations have 
mostly not been included so far in theoretical models. 
However, there are some systems which have been re- 
cently considered in this field where such correlations 
may play an important role. In the first place, RSOC has 
been observed in oxide heterostructures [l3|- Important 
effects of a RSOC have been experimentally observed in 
other transition metal oxides, particularly the iridatcs 
where the relevant degrees of freedom are the 5d electrons 
on the Ir^"*" ions, and where strong correlation effects are 



expected ^18'|. It should be also taken into account that 
correlation effects become increasingly important in sys- 
tems with low spatial dimensions. Interesting crossovers 
have been reported when moving from two-dimensional 
(2D) systems to quantum wires ^19j. In this sense, spin- 
orbit effects, predominantly of the Rashba type have been 
studied in quantum wires 2^ 21 1 . The presence of RSOC 
in g raphene and carbon nanotubes has been pointed out 
[22] and in fact a theoretical study with up to two inter- 
acting electrons has been done on a nanotube quantum 



dot [23 1 



In this Letter, an attempt of studying the interplay be- 
tween electron correlations and RSOC will be performed 
by considering a Hubbard model with a Rashba SOC. 
During the last two decades it has been realized the 
difficulties of tackling electron correlations theoretically. 
In studies of RSOC, electron correlations have been in- 
cluded only locally regarding quantum dots 24, 25] or on 
effective onedimensional (ID) systems using Tomonaga- 
Luttinger liquid approach [26| . For this reason, in the 
present Letter a two-leg ladder, which may be the nar- 
rowest quasi- ID geometry on which a spin-Hall effect 
could take place will be considered. The Hubbard- 
Rashba model on a two-leg ladder is suitable for treat- 
ment using essentially exact, unbiased, computational 
techniques, which have been enormously valuable in the 
field of strongly correlated electron systems. Particularly 
for mesoscopic systems, the full quantum many-body in- 
terplay between different terms of the Hamiltonian can 
be reliably captured. Then, the main goal is to determine 
and understand the behavior of the spin accumulation as 
a function of the RSOC and the Hubbard repulsion U. 
The interplay between RSOC and U with respect to a 
possible metal-insulator transition [27| will not be dis- 
cussed in the present work. 
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FIG. 1: (Color online) (a) The two-leg ladder system here 
considered; (b) a two-leg ring pierced by a magnetic flux; (c) 
previously study 2D rings; (d) setup used in DMRG calcula- 
tions. 



Formulation. The Hamiltonian of the Hubbard- 
Rashba model in a square lattice is defined as: 



H 



</,m>,(T / 

I 

cly^^c,,^) + H.c] (1) 



where the notation is standard. The first two terms cor- 
responds to the usual Hubbard model and the last term 
contains the RSOC 



28|. The couplings <o and Vso are 



normalized in such a way that + V^q — 1 (which is 
henceforth taken as the unit of energy). The longitudi- 
nal (transversal) direction of the ladder corresponds to 
the cc-axis (y-axis), as shown in Fig. H^a). The RSOC 
in ([T]) corresponds to an effective Rashba electric field 
En along the z-axis, i.e. perpendicular to the plane of 
the ladder. Open boundary conditions (OBC) are ap- 
plied in the transversal direction. The total number of 
lattice sites is = 2 x L. Hamiltonian ^ corresponds 
to the isotropic ladder, in which couplings are the same 
along both directions. The present study is limited to the 
quarter filled system (n = 0.5). Weak-coupling analysis 
[29| has shown that the Hubbard model on the isotropic 
ladder at n = 0.5 is at the boundary between a metal 
(CISO) and an insulating (COSl) phase [s^]. Exactly at 
the boundary, where additional metallic phases are pre- 
dicted to be stable, both this study as well as numerical 
studies [sil [s^l and the present results point to a metal- 
lic phase, even for the values of [/ > here examined. 
As soon as Vso is turned on, the occupied bonding band 
will split leading to a four-point Fermi surface and the 
system will behave as a metal. In all cases, the length L, 
equal to the total number of electrons, was set equal to 
Am {m integer). 

In order to induce persistent currents, a magnetic fiux 
piercing the ring formed by a closed ladder has to be ap- 
plied (Fig.[IJb)). This magnetic fiux $ is included in the 
Hamiltonian by the usual Peierls factors, which in mo- 
mentum space is equivalent to replacing by fc^^ -I- $ /L 
(units in which e = c = h — 1 have been adopted). 
These Peierls factors are included in both the hopping 
and SO terms of The resulting Hamiltonian is stud- 
ied using exact diagonalization (ED) for U > 0. For the 



noninteracting case, (7 = 0, by working in momentum 
space, ladders with L up to 6400 have been considered 
[40|. Due to the RSOC the total is not conserved 
hence all possible values of total have to be included 
in the Hilbert space. In order to reduce the Hilbert space 
dimension, translation symmetry along the longitudinal 
direction has been implemented thanks to the use of peri- 
odic boundary conditions (PBC) along the x-axis. Notice 
that the ladder ring shown in Fig. [Ijb) is different from 
the 2D ring (Fig. [Ijc)) that has been intensively stud- 
ied in the non-interacting case [ssl - lsij . Studies including 
Coulomb interactions for these rings have been limited 
to few electrons [36| . 

The total current is defined as J{^) = '^^g^'^^ where 
Eo{^) is the ground state energy. The contribution from 
the hopping term in ((T|), the hopping or spin-conserving 
current, is then given by the usual current operator along 
the legs, jh = 'itoJ2i,cri'4+x,a^i,'^ ~ H.c). A similar cur- 
rent operator can be easily derived for the contribution 
from the RSOC term in ([T|), that is, the SO or spin- 
fiipping current, jso- Using the Feynman-Hellmann the- 
orem, the total current is given by LJ — jh + jso- The 
spin polarization due to the spin-Hall effect is defined as 
AS*^ =< S'f — 5*1 > /2, where Sj is the total z-component 
of the spin on leg j. All physical quantities are ground 
state averages, and hence, they are functions of the mag- 
netic flux $. In order to determine the effects of the 
Rashba and Hubbard couplings on these quantities, the 
reasonable choice of looking at their maximum values as 
a function of the flux $ was adopted. 

The second setup considered in this work is shown 
in Fig. [Ijd), which corresponds to OBC along the x- 
axis. Here, at a given time the two ends of a lad- 
der are connected to a voltage bias AV generating a 
time-dependent regime which is studied by DMRG. Al- 
though various sophisticated approaches have been pro - 
posed within DMRG to deal with this time evolution 37 - 
|39[, the simple "static" approach is implemented in the 
present work which is enough to reliably capture the main 
physical features in the small lattices studied. This is a 
non-equilibrium process, but a small value AV = 0.01 
was adopted. In this case, the total current is computed 
as J{t) — where Ni is the total charge in the left 
half of the ladder, and t is the time measured in units 
of the time increment At = 0.1. Jmax is computed as 
the average of J{t) between the first two peaks of its 
time evolution, which can be precisely determined with 
a relatively small number of retained states M within the 
static approach. The spin-conserving current jh{t) was 
computed at the central links of the ladder. jh,max was 
computed using the same criterion as for Jmax- 

Results. For the case of persistent currents. Fig. [T{b), 
results for the noninteracting [/ = case are shown in 
Fig. [2] For various lattice sizes ranging from L = 8 to 
6400, it can be observed the maximum values of the spin- 
fiipping current (Fig. [H^a)) and the spin accumulation 
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FIG. 2: (Color online) Maximum value of (a) spin-orbit 
current, and (b) spin accumulation AS^, as a function of 
Vso /to for the noninteracting case on various two- leg ladders, 
n = 0.5, tight-binding calculations. 
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FIG. 3: (Color online) Maximum value of (a) total current, 
(b) hopping current, (c) spin-orbit current, and (d) spin ac- 
cumulation AS^ as a function of Vso /to for various values of 
U. Results for the 2x8 ladder, n = 0.5, obtained with ED. 



(Fig. [2jb)) occur at Vso/to ~ 0.5. The curves are quite 
noisy due to various level crossings as a function of the 
flux $ but nonetheless the ratio jso,max/ ^S^^^ follows 
a smooth monotonous behavior as a function of Vso /to 
for large L, as shown in [40| . Taking into account this 
noisy behavior, AS^^^ is quite independent of L. 




FIG. 4: (Color online) Maximum value of (a) total, (b) hop- 
ping, and (c) spin-orbit currents, divided by AV, and (d) spin 
accumulation AS^ as a function of Vso/to for various values 
of U. Results for the 2 x 20 ladder, n = 0.5, obtained with 
DMRG. 

Let us discuss the correlation effects for the case of 
persistent currents, for < [/ < 8. Fig. [3] shows the re- 
sults for the maximum values of the total current, J, the 
hopping and spin-orbit contributions to the current, jh 
and jso I and the spin-leg polarization AS*^ on the 2x8 
ladder at n = 0.5. As it can be seen in Fig. [2Ja), the 
total current for a fixed value of Vso /h decreases with 
increasing U. This result is expected for Vso = 0, i-e. 
for the pure Hubbard model, where it is well known that 
the conductance G decreases with increasing [/, as in the 
one-dimensional Hubbard model presenting a Luttinger 
liquid behavior where G ~ Kp |41| , as examined for lad- 



ders in [42|. The maximum value of the hopping current, 
shown in Fig. |3l^b), follows roughly the same behavior 
with U and Vso/to as the total current since it is much 
larger than the spin-flipping contribution. The dip in J 
and jh at Vso/to ~ 0.25 (Fig.|3l[a),(b)) is just afinite size 
effect. On the other hand, the behavior of the maximum 
value of the SO current is quite interesting (Fig. [3l[c)). 
For all values of U, jso, max increases with Vso/to un- 
til reaching a maximum which depends on U , and then 
starts to decrease becoming eventually negative, indicat- 
ing that its direction is opposed to the hopping current, 
by convention. The most important results correspond 
to the maximum value of the spin-leg polarization shown 
in Fig. [Sjd). One should first notice that AS^^^ in- 
creases as a function of the Coulomb repulsion U for all 
values of Vso/to- In addition, its maximum value near 
yso/to — 0.5 observed for the noninteracting case is pre- 
served for finite U . The dependence of these properties 
as a function of the magnetic flux, for Vso/to = 0.5 and 
1 , is shown in [iO] ■ 

Results for the voltage bias setup (Fig. [TJd)), for the 
2 X 20 ladder, are shown in Fig. SI The behavior of the 
maximum values of the total, spin-conserving and spin- 
flipping currents, and of the spin accumulation as a func- 
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FIG. 5: (Color online) (a) (S^)^ as a function of Vso/to for 
[/ = 0, 2, 4, 6, 8 from bottom to top for the 2 x 8 ladder PBC at 
zero flux (circles) and for the 2 x 20 ladder OBC (diamonds), 
(b) Peak value of the magnetic structure factor ^^^(q) as a 
function of Vso/to for U = 0, 2, 4, 6, 8 from bottom to top for 
the 2x8 ladder PBC at zero flux (circles, stars) and for the 
2x20 ladder OBC (diamonds, crosses) if q = (vr, tt) (other) (c) 
Maximum value of the spin accumulation AS^ as a function 
of (5^)^ at zero flux for various values of Vso/to indicated in 
the plot, for the 2 x 8 ladder PBC. 



tion of Vso/to for various values of U, are qualitatively 
strikingly similar to those shown for the setup involving 
persistent currents (Fig. [3]) in spite of the quite different 
background physics involved. Notice that actually the 
maximum values of the currents divided by the voltage 
bias Ay are shown, which correspond to the conductance 
G up to a factor of 2tt. Comparison with results for the 
2 X 16 do not show appreciable finite-size effects. Details 
of the calculations are provided in . 

Some possible explanations of the behaviors shown in 
Figs. [3] and [4] can be obtained by examining the depen- 
dence of various magnetic properties with the parameters 
U and Vso/to- In Fig. [5]Ja) the average value of (5"^)^, 
a measure of the magnetic moment per site [40| . as a 
function of Vso/to and for various values of U is shown. 
As expected, (S'^)^ increases with U, as it is well-known 
for Hubbard-like models, and it is almost independent of 
Vso/to- Results for the 2x8 ladder PBC at zero flux 
are virtually indistinguishable with the ones obtained for 
the 2 X 20 ladder OBC using DMRG. More interesting 
are the results for the peak value of the static magnetic 
structure factor S'^^(q), which are shown in Fig.lSJb) as a 
function of Vso/to and for various values of U. For large 
values of Vso/to, S^^{q) has a peak at q = (tt, tt), while 
for small values of Vso/to, S^^{q) is maximal at other 
values of q, all of them with tt, and in most cases at 
q = (0, tt). In the first place, as in Fig.ISja), it is expected 
that S'^^(q) is enhanced by U. What is unexpected is 
the strong dependence of S'^^(q) with Vso/to- S^^{TT,n) 
increases with increasing Vso/to, and ^^^(q ^ (tt, tt)) in- 
creases with decreasing Vso /to- This is a combined effect 



of Vso/to and U since for U = 0, S'^^(q) is constant and 
clearly the slope of ^^^(q) with Vso/to increases with U. 
Results for the 2 x 20 ladder show a stronger dependence 
than the ones corresponding to the 2x8 ladder and this 
could be explained by the fact that in DMRG S'^^(q) cal- 
culations, the only included spin-spin correlations are the 
ones measured from one of the central sites of the ladder. 

The increase of the z-component of the magnetic mo- 
ment can at least partially explain the larger values of the 
spin accumulation observed as a function of the Hubbard 
repulsion for a fixed value of Vso / to - This is shown in 
Fig. [Sfc) where it can be observed an almost linear de- 
pendence of AS'^j,^ with (S^)^ except for a slight upward 
curvature for Vso/to = 0.5 and a slight downward cur- 
vature for Vso/to ~ 0.75. The peak of the magnetic 
structure factor at q = (7r,7r), indicating the presence 
of an antiferromagnetic (AFM) order could in turn pro- 
vide some insight on the suppression of the spin-flipping 
current jso for large Vso/to and U, shown in Figs.[3][c) 
andSfc). For an AFM order, two electrons occupying 
nearest neighbors sites along one leg of the ladder would 
likely have spins with opposite projections on the z-axis. 
For this configuration, a spin-flipping hopping would be 
impossible while there would be no problem for a spin- 
conserving hopping (except for a cost U). Now, once 
a io-hopping has occurred, creating a doubly-occupied 
site and leaving behind an empty site, a Vso-hopping is 
now possible from the doubly-occupied site to the empty 
one. Hence the spin-flipping hopping takes place in a di- 
rection opposite to the spin-conserving one thus possibly 
explaining the observed change of sign of the spin- flipping 
current. 

Conclusions. To summarize, we have considered per- 
haps the simplest system where some features associ- 
ated to Rashba spin-orbit coupling can be observed [3] . 
This two-leg ladder system is suitable to be studied by 
unbiased, essentially exact, computational techniques, 
which allow to consider electron correlations on an equal 
foot with spin-conserving and spin-flipping hoppings. At 
quarter-filling it was observed that spin accumulation is 
maximal at an intermediate value of Vso/to, indepen- 
dently of U. By increasing a repulsive U, the spin-flipping 
current is suppressed and it becomes negative for large 
Vso and U, while spin accumulation is enhanced. These 
two last features can be explained by the increase of (S^)^ 
and AFM correlations, the latter being a combined ef- 
fect of Vso and U. The present approach could be ex- 
tended to wider strips which could in principle be studied 
by DMRG, and also by variational and diffusion Monte 
Carlo (plus a fixed-node approximation), and to differ- 
ent electron fillings. In addition, some subtle and so far 
controversial properties, such as spin currents, could be 
studied numerically on this class of strips including elec- 
tron correlations. Although the present study does not 
include any dissipative mechanism, an enhancement of 
spin accumulation by electron interactions together with 
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a reduction of the overall current are positive features 
that could be taken advantage in spintronic devices. 
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Supplementary material for: Spin 
polarization in the Hubbard model with 
Rashba spin-orbit coupling on a ladder 

Jose A. Riera 

Institute de Fi'sica Rosario y Departamento de Fisica, 
Universidad Nacional de Rosario, Rosario, Argentina 

The tight-binding part of the Hamihonian (1) in the 
main text is, in momentum space: 
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in the basis: (ci,k,ti ci,k,_|„C2,k,t7 C2,k4,): where "1", and 
"2" are the two legs of the ladder. 
For the general spatially anisotropic ladder: 

A = ^2t|| cos 

B = -2jV^5o,|| sinfc^ (S2) 

In the present case, t\\ = t j_ = to, Vso,\\ = Vso,!- = Vso- 
In the tight-binding formalism, the currents J, jh , and 
jso are independently computed for each value of the 
magnetic flux and hence the relation LJ — jh + jso 
provides an internal check of the calculations. Similarly, 
in the ED and DMRG calculations, J and jh are indepen- 
dently computed, hence for Vso — the relation LJ = jh 
has been checked. 




FIG. SI: Ratio of the maximum values of spin flipping cur- 
rent to the spin accumulation as a function of Vso /to for the 
noninteracting case. Persistent currents setup. Sizes of the 
ladder rings are indicated in the plot. 

Fig. IS II was built from data of Fig. 1 of the main text. 
Within numerical accuracy, there is a local minimum at 
Vso/h = 1 with the value jso/^S'' w 0.848. 
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FIG. S2: (a) Ground state energy relative to the value at zero 
flux, (b) spin conserving current, (c) spin flipping current, and 
(d) spin accumulation, at Vso /to = 0.5, on the L = 8 ladder 
for various values of U , and for the L — 6400 ladder for the 
noninteracting case, as a function of the magnetic flux. The 
values of E{^) — E{0) for L — 6400 have been multiplied by 
400. 
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FIG. S3: Same as Fig. [S2] but for Vso Ac = 1- Notice the 
inversion of the spin-flipping current Jso for (7 = 4. 



In Figs. [S2] and IS3I the evolution of various physical 
properties with the magnetic flux is shown for Vso /to — 
0.5 and Vso /to — 1-0 respectively for the L ~ 8 ladder 
and various values of U. Fig. lS2l shows an almost negligi- 
ble difference with the results obtained for the noninter- 
acting case with the L = 6400 ladder. On the other hand, 
in Fig. [S3] very large differences can be observed. What 
is remarkable in this plot is that the spin-flipping current 
jso (Fig- [S3l c)) is inverted for [/ = 4 thus becoming 
opposite to the spin-conserving current jh (Fig. ISSfb')'). 

In Fig. [S4l some typical time evolution of the total cur- 
rent and of the spin accumulation are shown for various 
values of U and Vso /to- It is well-known that the time 
evolution of the current follows an oscillatory behavior 
due to the interplay between kinetic and potential en- 
ergies. As the number of retained states increases, or 
by adopting an adaptive scheme, the current reaches a 
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FIG. S4: Time evolution of the total current J divided by AV 
(upper panels) and of the spin accumulation (lower panels) for 
(7 = (left panels), [7 = 4 (central panels) and U = 8 (right 
panels), and various values of Vso/to, on the 2 x 20 ladder, 
obtained by DMRG. Only the first half-period of the time 
evolution is shown, for the next half-period both J and AS" 
reverse their sign. 




FIG. S5: (a) Weight of each sector in the ground state, 
and (b) contribution of each 5*^ sector to the total average 
(S^)^ for U = 0. (c), (d), same for U ^ 8. 



plateau during each half-period of the time evolution. 
The value of the current at the plateau divided by the 
voltage bias AV is equal G/27r, where G is the con- 
ductance (remember that units where e = h = I were 
adopted). In the present case, the adopted criteria of 
measuring Jmax as the average of J{t) between the first 
two peaks leads to a 7% overestimation of G with respect 
to its exact value at ?7 = Vso = 0. 



(S*^)^ and the static magnetic structure factor are de- 
fined as usual from the z — z spin correlations: 
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With this definition, is easy to find that, for the non- 
interacting case, (S^)"^ = nf{l — nj,) -I- (t^i), where n-f 
(rii) is the density of electrons with t {],) spin in a state 
with fixed total 5^ Note that for 5^ = 0, (5^)^ = 0.375 
which is the result shown in Fig. 5(a) of the main text, 
in spite of the fact that is not a quantum number 
for the model here considered, and the ground state is 
a combination of subspaces with all possible S^. Also, 
(5^)^ n = 0.5 when [/ — !> oo since doubly-occupied 
sites become forbidden. Fig. IS5I shows the weight of each 
sector in the ground state, and the contribution of 
each sector to the total average (S^)'^. 



